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Abstract 

We consider a hydrogen atom confined in a thick brane embedded in a higher-dimensional space. 
Due to effects of the extra dimensions, the gravitational potential is amplified in distances smaller 
than the size of the supplementary space, in comparison with the Newtonian potential. Studying 
the influence of the gravitational interaction modified by the extra dimensions on the energy levels 
of the hydrogen atom, we find independent constraints for the higher-dimensional Planck mass in 
terms of the thickness of the brane by using accurate measurements of atomic transition frequencies. 
The constraints are very stringent for narrow branes. 

PACS numbers: 04.50.-h, 11.25.Mj, 04.80.Cc, 32.30.-r 
Keywords: extra dimensions, branes, atomic spectroscopy 


1 


I. INTRODUCTION 


Since the beginning of the past century, theories of higher-dimensional spacetime 
have attracted attention of physicists mainly motivated by the Kaluza-Klein theory, which 
postulated the existence of a compact extra dimension with the size of the Planck length in 
order to formulate a unihcation theory of electromagnetism and gravitation. More recently, 
braneworld models, originally conceived as a framework to solve the hierarchy problem, 
have vigorously renewed the interest in this subject ij-l^. In the braneworld scenario, our 
observable universe is a submanifold isometrically embedded in a higher-dimensional space 
and, unlike the Kaluza-Klein model, the size of the extra dimensions can be much larger 


than Planck length or even inhnite jd] without introducing any irremediable conflict with 
experimental data. 

The main feature of the braneworld scenario is that the standard model helds are conhned 
in a submanifold (3-brane), if they are not excited above a certain energy scale, which should 


be greater than TeV order j^. On the other hand gravity has access to extra dimensions. 
Despite this freedom, due to the existence of a bound zero-mode, the gravitational force re¬ 
covers its four-dimensional behavior for distances much larger than the characteristic length 
scale, £, of the supplementary space. In models with compact dimensions, £ corresponds to 
the size of the extra dimension, while in models of non-compactihed extra dimension, i is 
related to the curvature radius of the ambient space. Below this length scale, the influence 
of extra dimensions on the gravitational force become signihcant. It happens that i can be 
much greater than the scale, L^, in which the conhned helds feel the ehects of the extra 
dimension directly. In the original ADD model, < 10“^® m [l|, while, for branes inspired 
in string theories, is much smaller j^. 

In sight of this, the search for traces of extra dimensions in the gravitational interac¬ 
tion arises as a paramount issue which can be empirically investigated. Direct tests of the 
gravitational force in microscopic domains, based on torsion balance experiments, found no 
deviation from the inverse square law, constraining, in this way, the radius R (= £/27r) of the 


Q. 


extra dimension. Considering the original ADD model [li , it is founc 
is the tightest bound for models with only one extra dimension 
mensions, the most stringent constraints come from astrophysics 
particle collisions EQ. 


hat i? < 44 pm, which 
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Here we study the effects of the extra dimensions in a bound system, in order to find 
independent constraints for the number and the size of the supplementary space. It is 
reasonable to expect that the extra dimensions modify significantly the behavior of gravity 
in distances shorter than £. If the gravitational field in the weak-field limit obeys Gauss law, 
then the gravitational potential of a point-like mass m has the following behavior: 


V(r) = 


— for r S> ^ 
for r < f 


( 1 ) 


where Gn is the gravitational constant in the higher-dimensional space with n extra di¬ 
mensions. The exact relation between Gn and the Newtonian constant G depends on the 
extra-dimensional volume, which, on its turn, depends on the topology and curvature of the 
supplementary space. However, in magnitude order, it is expected that Gn ~ Gt^. 

In atomic systems the usual gravitational interaction is negligible. The average potential 
gravitational energy of a hydrogen atom in the ground state is of the order of Gmpirte /~ 
10-38 where Oq is the Bohr radius. However, due to effects of the extra dimensions, the 
potential energy increases by a factor of the order of for short distances (r < £). Thus, 

depending on the values of i and n, the gravitational energy could in principle modify the 
energy spectrum of the hydrogen significantly. Therefore, the highly accurate measurements 
of atomic transitions can be used to put empirical constraints on models of extra dimensions. 

Treating the gravitational potential as a perturbation of the Hamiltonian of the hydrogen 
atom, we find that the energy shift caused by extra dimensions in a certain quantum state 
ip is proportional to iPg average of r“G+i) that state. It happens that, for 

n > 2, this average diverges for S'-levels. This problem is connected with the fact that the 
short-distance behavior of the brane-to-brane graviton propagator is not computable in the 
original ADD model, even in the tree level, suggesting that it may depend on the ultraviolet 
details of a fundamental quantum theory of gravitation, as it was pointed out in Ref j^, 
where the effects of extra dimensions on high-energy particle collision via the gravitational 
interaction were investigated. 

The first attempts to constrain extra-dimensional parameters by using atomic spec¬ 
troscopy have tried to circumvent the divergence problem by establishing a cut-off radius 


14l-ll8j. However, proceeding in that way, the constraints become dependent on this arbi¬ 


trary parameter. Another way to avoid the divergence problem is by considering the effects 


3 






of the modified potential on molecular spectroscopy instead of atomic spectroscopy |19 |. 
However, this approach gives much weaker constraints. 

In the classical level, the divergence problem arises because the matter distribution con- 
hned in the brane is singular from the perspective of the ambient space. It happens that 
such delta-like conhnement is just an idealization that hides the internal structure of the 


brane. The relevance of the brane substructure to regu 


fective theory has been observed and studied before 


20 


ate divergences of a thin brane ef- 


21| . In Ref 20[, for instance, these 


divergences are treated by introducing a renormalization prescription. Here we follow a dif¬ 
ferent procedure which relies on the fact that, in the thick brane scenario, the wave function 
of localized particles has a certain width a in the extra-dimensional space. As we shall see, 
when the width is taken into account, the average is hnite and depends on a. As 

the value of width is bounded by the thickness of the brane, then our analysis provides a 
joint constraint of the higher-dimensional Planck mass and brane thickness. 


II. AN ATOM IN A THICK BRANE 


In the context of string theories, the brane may be considered as inhnitely thin. However, 
in the framework of the held theories, it may have a thickness. For instance, if we consider an 
ambient space of (4-|-^dimensions, a 3-brane may be represented as a domain wall solution 
of some scalar held |22|. A Yukawa-type interaction between Dirac spinor helds and the 
scalar held can localize matter in the core of the domain wall. Due to this interaction, the 
state of zero-mode of the spinor held is given by the following wave function: 


T (x, z) = exp 


-h / (t>o{y)dy 
Jo 


V’(x), 


( 2 ) 


where h is the coupling constant, ■0 (x) represents a free spinor in the (3 -|- l)-dimension, 
00 = Tj ta.Vih {z/e) is the scalar in a domain wall conhguration interpolating between two 
vacua ±7] of the scalar held. With respect to the transversal direction z, the wave function 
has a peak at 2 ; = 0 and exponentially falls as we move away from the center of the domain 


wall. The quantity e: dehnes the thickness of t 


le brane and must be smaller than 10 


, 3 - 


m, 


according to current experimental constraints 

Following a similar reasoning, it is possible to devise a conhnement mechanism for matter 
in topological defects of greater codimensions. By generalizing the above result for a space 


4 


















with n extra dimensions, we can consider that, from the phenomenological point of view, 
the wave function of localized particles can be written as 'h (x, z) = x(z)'0(x). For the 
sake of simplicity, we will assume that the extra-dimensional part y is a Gaussian function 
around the center of the brane with a standard deviation a: 

n/4 

2 ] 


X{z) = 




n n 

■ zf 




( 3 ) 


2 = 1 


Of course the width of the wave function in the transverse directions should satisfy a < e < i. 
At this point it is important to emphasize that, although we have chosen a Gaussian prohle 
in the transverse direction, x can be any normalizable function of z, as we shall see later. 

Now we can describe the gravitational interaction between the proton and the electron in 
this scenario. In the linear regime of gravity, there exists a coordinate system (a gauge) in 


which the gravitational potential satishes the Laplace equation in vacuum 


23 |. Thus, in the 


hrst order of the exact potential produce by a point-like mass in a 


ligher-dimensional 


space can be obtained, at least formally, for any compact topology For instance, if 


the supplementary space has a topology of a n—dimensional torus with size £, then the 
gravitational potential of a mass m at the point R = (x, z) of the ambient space is given by 

Grim Grim 


G R = 


^n+l 




R-R' 


/ in+l ’ 


( 4 ) 


where R' = £ {0,0,0 ki...,kn) and each ki is an integer number. The potential (H|) is a 
solution of the (3 -I- n)-dimensional Laplace equation with the appropriate boundary con¬ 
ditions. Indeed, the form of vector R,; ensures that V (R) is periodic with respect to the 
extra-dimensional directions as required by the topology of the supplementary space. It is 
interesting to note that, from the perspective of the space V (R) can be viewed as a 

superposition of potentials generated by the real mass m at the origin and by copies (mirrors 
images) of the mass localized by the vectors R'. 

In the far zone, i.e., for |x| >> £, (|1]) reproduces the Newtonian potential. Moreover 
the corrections due to the extra dimensions can be written as 1/r (l -|- ae [2^. This 
Yukawa-law form is vastly employed to constrain parameters of higher-dimensional models, 
by using data from the torsion balance experiments 1]|. On the other hand, in a bound 
system, such as the atom, the major contribution comes from the hrst term of (jlj), at least 
if the atom is at the lowest states. The hrst corrections come from the mirror images closest 
to m. The potential of each one of these 2n hrst neighbors is less than Gnm/£‘^~^^ ~ Gmj£, 
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which is negligible in comparison to the first term, as we can verify by calculating this term 
with the extra-dimensional estimated size, i, that we obtain in Figure 2. For a detailed 
estimate of the contribution of the mirror images see, for instance, the appendix of Ref. 
25|. 

Thus, for our purpose, we can approximate the gravitational potential of a point-like 
source by the function which is proportional to the Green function of the 

Laplace operator in the flat higher-dimensional space Therefore, in this approxima¬ 

tion, the potential produced by the proton, assuming that its mass nip is distributed on the 
spatial extension of the nucleus, is given by 

= ( 6 ) 

where the mass density is p = nip and Tp (x, z) = y (z) ifjp (x) is the higher-dimensional 
wave function of the proton. We assume that the mass is uniformly distributed inside the 
nucleus, so the 3-dimensional part 'i/'p (x) is constant in the spatial extension of the nucleus 
and zero outside. 

In the hydrogen the average kinetic energy of the electron is a small fraction of its rest 
energy. Thus, in a first approach, we are going to consider the non-relativistic picture of the 
atomic system. In the classical framework, the gravitational interaction between electron 
and proton is governed by the Hamiltonian Hq = nigcj). Now assuming that Hq is a small 
term of the hydrogen Hamiltonian, it follows, from the perturbation method, that the energy 
shift corresponding to a certain atomic state, in the first order, is given by 

6E^ = |Te|'me0(R)d3+”R, (6) 

where Te (x, z) is the higher-dimensional wave function of the electron (more precisely, the 
reduced particle) which comprises the extra-dimensional part y (z) and the usual solutions 
'ipe (x) of the Schrodinger equation for the hydrogen atom. 

Therefore, the extra dimensions modify the energy spectrum of the hydrogen by means of 
the gravitational interaction between electron and proton. Considering the wave functions of 
IF and 2S states {i/jis = 1/ \/ Traj] exp {—r/ao) and ‘ijj 2 S = 1/ v^STrag (1 — r/2ao) exp (—r/2ao), 
respectively), the shift in the frequency of the 2S — IS transition, Az/g = \ 6 E 2 s — SEis\ /h, 
can be calculated from ([6]) and compared with the experimental value. 
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III. RESULTS AND DISCUSSION 


The empirical value of the 2S — IS* transition frequency in the hydrogen atom is 


/exp = 2466061413187035 Hz with an experimental error 5exp = 10 Hz 


261. The theoretical 


prediction (based on well-established four-dimensional physics) agrees with the experimental 
value up to the order of the theoretical error (Sth), which corresponds to 32 kHz and is related 
to uncertainties in measurements of the proton radius 27|. Thus, to be consistent, any effect 


of extra dimensions should be lesser than the combined errors. The constraints on the pa¬ 
rameters of the extra dimensional theory are obtained by imposing that Auq < \J^th + ^exp- 
The numerical analysis of this condition is summarized in Figure 1. It establishes a lower 
bound for the fundamental Planck mass of the higher dimensional space Mjy in terms of a. 


According to Refs. 


23j, the fundamental Planck mass is related to by the formula 




= Qn{fi/c)^ hc/Gn, where = (27r)’^ P /[{n + 2)27r(”+^)/^]. In Figure 1, the 
regions below the curves are excluded. It is worth mentioning that the width of the particle 
wave function in the transverse directions, a, is smaller than the thickness of the brane (e). 
Therefore, the analysis shows that the constraints are tighter for thinner branes. 

The cases n = 1 and n = 2 are not shown because the corresponding constraints are 
very weak. For n > 3, the major contribution of the gravitational energy comes from the 


-19 


m. 


integral in the interior of the nucleus, when we consider realistic branes with £ < 10 
The leading term is of the order of GnTUpme/. This dependence on a explains why 
narrow branes provide very stringent limits on the fundamental Planck mass. 
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Figure 1. Lower bound for the higher-dimensional Planck mass M/j (expressed in natural units) 
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in terms of (T(the wave function width in the transverse directions).The areas below the lines are 

excluded. 


In the original ADD model, the snpplementary space is a flat n—torus with radius R. 
The relation between R and the higher-dimensional Planck mass Mq (given by G~^ = 
whereAT^"" = he], see jb]) can be used to constraint the radius 

of the extra dimension. Figure 2 shows upper limits on the extra dimension radius for 
n = 3,4, 5 and 6. The regions above the curves are excluded. 
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Figure 2. Upper limit for the radius of the extra dimensions in terms of afor n=3, 4, 5 and 6. 

The regions above the lines are excluded. 


As mentioned before, we have admitted that the wave function in the supplementary 
space has a Gaussian prohle ([3]) for the sake of simplicity. Nevertheless, it is important to 
emphasize that this assumption does not play an essential role in our analysis. Actually, 
our results do not change, even if we consider other prohles, provided that the parameter a 
is dehned by: 


1 _ r(n/2) 
^ = r (M) 


\Xp{Zi)f \Xe {Z2)f 


d"'Zid"'Z2, 


(7) 


I ^ ^ \m 

Id - d| 

where m is a positive integer that should satisfy the condition m < (n — 1) and T stands for 
gamma function. When the prohle is Gaussian, this parameter a corresponds to the width 
of the Gaussian function. 

At this point, we should stress that the above constraints were derived based on the clas- 

r 

sical regime of the gravitational interaction. Nevertheless, as it was emphasized in Ref. [6[, 
quantum-gravity effects may become relevant in a length scale of the order of Ijj (the higher 










dimensional Planck length, given hy Id = i^/c) ) or even in a greater scale, depending 

on the quantum-gravity theory, not yet known. If this is the case, then unexpected effects 
may emerge and could distort or even overshadow the classical results we have investigated 
here. 

n the lack of a quantum-gravity theory, let us explore some possibilities. According 
, quantum corrections to the gravitational potential energy may be estimated by 


to 


treating General Relativity theory as an effective theory. The classical potential energy, in 
three-dimensional space, is given by the Newtonian term GMm/d, where d is the distance 
between particles with mass M and m. The quantum contributions are smaller by a factor 
of the order of [Ip/df', where Ip is the Planck length of the ordinary three-dimensional 


28l |. In the higher-dimensional case, the classical potential is and the 


space 

quantum corrections would be of the order of {In/d)^'^'^, based on dimensional analysis. 
Now, considering the hydrogen atom, it is instructive to define an effective extra-dimensional 
distance, dg//, between proton and electron, though they cannot be considered as point-like 
particles in this bound system. By writing the atomic gravitational energy as G^mpme/d”^^, 
it follows that d^^j ~ cT"“^ag in the ground state. Notice that dg// is a kind of geometric 
average of the wave packet widths in the transversal and parallel direction of the brane. Now, 
from the lower constraints on Md, we can estimate upper bounds for Id as a. function of a. 
By a direct calculation, we can verify that the ratio {Id/ depends on n and cr, but, for 
any dimension and for any value of a investigated here, it is smaller than 10“®. Thus, if dg// 
is the relevant characteristic length scale that regulates the gravity behavior in this system, 
then we may expect that the classical contribution will be the leading gravitational influence 
in this context. However, as the fundamental quantum-gravity theory is not known, only 
experiments can decide on this question. 


IV. COMPARISON WITH OTHER CONSTRAINTS 

Experimental bounds for the fundamental Planck mass are determined from many areas 
of Physics. As we have mentioned before, when the codimension is greater than 2, some 
of the most stringent constraints of Md are established by high-energy particles collisions. 
Data analysis from monojet events in proton-proton collisions at the LHC gives the fol¬ 
lowing lower bounds for Md in TeV/c^: 4.38 {n = 3), 3.86 (n = 4), 3.55 (n = 5) and 3.26 
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(n = 6) 29|, l30|. Examining Figure 1, we may say that, in this thick brane scenario, unless 
unpredictable quantum-gravity effects suppress the classical result, spectroscopy could give 
stronger constraints in the case the confinement parameter, a, is small. 

However, as a has an unknown value, it is interesting to investigate the inverse problem, 
i.e., to estimate the maximum infiuence of the gravitational energy on the hydrogen spec¬ 
trum, considering the current constraints on Mr, given by the LHC. In Figure 3, we show 
the gravitational contribution to the hydrogen (25 — 15) transition, AEg, as a function of 
the confinement parameter, a. 
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Figure 3. The energy gap between the 25 — 15 states in the hydrogen atom due to 
higher-dimensional gravitational interaction, assuming LHC constraints for M^t. The horizontal 
lines, labeled with and SEthSi^^e the current experimental precision and theoretical 

uncertainty, respectively, in the 25 — 15 transition. 


Here the gravitational energy AEq is confronted with the current theoretical uncertainty 
SEth and current experimental error 6 Eg^p. For each value of n considered, there is an 
interval in the cx-axis in which AEq is less than the theoretical uncertainty. However, even 
in this range, extra dimensions induce a huge amplification of the atomic gravitational energy 
in comparison with the ordinary three-dimensional case, which is of the order of 10“^® eV. 
Within these specific ranges of a. Figure 3 indicates how much the experimental precision 
should be improved in order to the hydrogen spectroscopy become capable of revealing traces 
of the supposed extra dimensions. For shorter a, the limit on AEq is still given by 5Eth- 
Once again, it is important to emphasize that this prediction is valid only if no quantum- 
gravity effects suppress the classical result. But, notice that, if quantum-gravity effects 
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reinforce the classical contribution, then signals of extra dimensions, though unpredictable, 
would be viewed earlier. 

New forms of constraints can be obtained by considering relations between a and Id in 
advance. Let us admit that Id is smaller than a by a factor x, i.e.. Id = xa, with x < 1. 
In this case, the condition that the energy shift between 2S and IS" states, due to the 
gravitational interaction, be smaller than the uncertainty 6 E yields the following constraint: 


c^Md > 


/ j- \ 3 2 2 

f nc\ rUpC rUeC 


I 

V^O 


1 1/4 


6 E 


X 


(n-2)/4 


( 8 ) 


where = 7'yn^n/8TT and coefficients 7 „, have the following values: 73 = 27r^/^,74 = 
47 r/ 3,75 = 7 r ^/^/3 and 70 = dvr/lb. Thus, for instance, if a = IOId, then the lower bounds 
for Md, assuming 6 E of about 10“^° eV, would be of GeV order, therefore, 10^ weaker than 
the constraints from LHC. Equation ([ 8 ]) also shows a weak sensitivity of Md with respect 
to the experimental precision. Indeed, if 6 E was reduced by a factor of 10^, then, Md 
would be improved 10 times, only. It is also interesting to compare these constraints with 
indirect limits obtained in lepton colliders, which seem to have a similar origin. Ref. jo] 
determined the maximum Md sensitivity which can be reached by studying the hnal state 
with photons and missing energy at an electron-positron collider. Considering ^/s = 1 TeV 
and integrated luminosity C = 200 fb“\ the predicted bounds for Mdc^ in TeV are: 4.0 
(n = 3), 3.0 (n = 4), 2.4 (n = 5), when the beam polarization is 90%. These limits are much 
greater than those from equation ([ 8 ]). 

Another way to hnd new constraints is, instead of hxing a direct relation between Id and 
cr, to require that the proton energy density in the extra dimensional space, p, be smaller 
than the critical density, pd = Mdc^, by some factor y. Admitting a Gaussian prohle 
in the transverse direction, we can estimate the proton density, ~ m^(? Rl, in the center 
of the thick brane. From the condition p /pd = y, we can write a in terms of y. By using 
this relation, we hnd constraints for the fundamental Planck mass in terms of y: 




(n-2)/8 


(9) 


-| (n-2)/8 


where /3„ = 0 ;”^* (4vr/3) ( 7 r/ 2 )"'^^ . Thus, taking y = 0.1, to make some estimation, 

we hnd Md > 1-6 GeV/c^ (n = 3) to Md > 76 GeV/c^ (n = 6 ), considering the current 
theoretical uncertainty dE^h. In comparison with the previous case, the constraints are much 
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more sensitive with respect the precision 6E. Indeed, if the condition ([9]) is calcnlated with 
the cnrrent experimental error, (5-Eexp, the constraints will be at least times stronger. 

Specihcally, for n = 6, we wonld have M^, > 32 TeV/c^, which is more stringent than the 
collider constraints. 


V. FINAL REMARKS 


In the thick brane scenario, it is expected that the matter and the standard model 
helds do not feel directly the effects of the extra dimension in a length scale less than the 
thickness of the brane e, while the modihcation in the gravitational held may arise in a 
scale i mnch greater than £. Taking this into acconnt, we have seen that the corrections 
of the gravitational potential dne to extra dimensions may give signihcant contribntion in 
a bonnd system. Tims, by nsing precision measnrements of 25* — 1 S' transition freqnency 
of the hydrogen, we have obtained new constraints for the higher-dimensional Planck mass 
in terms of a conhnement parameter, a, of the matter inside the brane, which shonld be 
smaller than the brane thickness. 

The constraints we hnd here are stronger for thinner branes. In fact, comparing onr 
resnlts with the cnrrent constraints from LHC data, which go from M^, > 4.38 TeV/c^ 


(n = 3) to Mo > 3.26 TeV/c^ (n = 6) 


29 




we may conclnde that, for some narrow 


branes, atomic spectroscopy conld impose very stringent constraints for the fnndamental 
Planck mass in higher-dimensional spaces. 

Finally, it is important to emphasize that, as the higher-dimensional Planck mass can be 
mnch smaller than the fonr-dimensional Planck scale, it is expected that qnantnm-gravity 
effects may arise mnch sooner in comparison with the traditional pictnre withont extra 
dimensions jo]. However, as the fnndamental theory is not known, the snpposed qnantnm 
effects are nnpredictable. On the other hand, the constraints from the hydrogen spectroscopy 
we hnd here relies on the classical behavior of gravity, thns, it is important to highlight that 
the present resnlts are validity only if the classical contribntions are not snppressed by 
qnantnm-gravity effects. 
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